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ABSTRACT 


In  recent  years,  spacecraft  have  become  increasingly  flexible.  The  design 
requirements  for  the  Bifocal  Relay  Mirror  spacecraft  include  controlling  jitter  at 
the  nanoradian  level.  Typically,  tight  pointing  requirements  require  high 
structural  stiffness,  at  the  cost  of  increasing  the  on-orbit  mass.  To  accomplish 
this,  while  minimizing  the  mass  of  the  spacecraft,  the  structure  will  have  some 
inherent  flexibility.  These  flexible  modes  will  interact  with  the  pointing  control, 
hence  affecting  the  payload  performance.  The  compensator  design  conducted  in 
this  thesis  achieves  order  of  magnitude  improvements  in  controlling  the  rate 
error,  hence  jitter.  This  thesis  starts  with  a  rigid  body  dynamic  model,  and 
develops  a  flexible  body  dynamic  model.  Once  the  model  is  developed,  the 
structure-controls  interaction  is  discussed.  Finally,  compensators  are  applied  to 
the  rigid  body  controller  to  mitigate  the  performance  losses  present  in  the  flexible 
body  system.  Through  classical  second-order  compensators,  the  angular  rate 
error  was  decreased  by  a  factor  of  ten.  Nonminimum  phase  notch  filters  and 
phase  lag  filters  were  used.  Ultimately,  the  phase  lag  filters  provided  the  best 
performance. 
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I.  INTRODUCTION 


Recently,  spacecraft  designs  have  become  more  flexible  in  comparison  to 
the  legacy  satellites,  which  are  more  rigid.  There  are  also  designs  that  involve 
multibody  motion,  such  as  a  moving  payload  on  a  slewing  central  bus.  The 
Bifocal  Relay  Mirror  combines  these  situations  for  the  purpose  of  relaying  a  laser 
from  a  ground  or  air  source  to  space  relay,  to  a  ground  or  air  target  or  receiver. 
The  target  may  be  cooperative,  in  the  case  of  optical  communications  or  other 
non-force  application  purposes,  or  uncooperative  in  the  case  of  high-energy 
applications. 

Increasing  the  precision  of  the  acquisition,  tracking  and  pointing  sequence 
demands  that  multibody  dynamics  and  flexible  dynamic  effects  be  combined  and 
addressed  from  a  control  design  perspective.  This  thesis  presents  the  rigid  body 
and  flexible  body  dynamics,  and  compares  the  uncompensated  and 
compensated  performance  of  the  spacecraft  using  a  flexible  plant.  Classical 
compensator  designs  are  reviewed  and  applied. 

The  design  requirements  for  the  Bifocal  Relay  Mirror  spacecraft  include 
controlling  jitter  at  the  nanoradian  level.  Typically,  tight  pointing  requirements 
require  high  structural  stiffness,  at  the  cost  of  increasing  the  on-orbit  mass.  To 
accomplish  this,  while  minimizing  the  mass  of  the  spacecraft,  the  structure  will 
have  some  inherent  flexibility.  These  flexible  modes  will  interact  with  the  pointing 
control,  hence  affecting  the  payload  performance.  The  compensator  design 
conducted  in  this  thesis  achieves  order  of  magnitude  improvements  in  controlling 
the  rate  error,  hence  jitter. 

Through  classical  second-order  compensators,  the  angular  rate  error  was 
decreased  by  a  factor  of  ten.  Nonminimum  phase  notch  filters  and  phase  lag 
filters  were  used.  Ultimately,  the  phase  lag  filters  yielded  the  best  performance. 
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II.  BIFOCAL  RELAY  MIRROR  RIGID  BODY  EQUATIONS  OF 

MOTION 


This  section  presents  the  derivations  from  Appendix  B  of  the  Bifocal  Relay 
Mirror  Technology  Development  paper  submitted  by  the  Naval  Postgraduate 
School  in  February  2002  for  the  National  Reconnaissance  Office  Directors 
Innovation  Initiative  competition.  It  provides  the  derivations  of  the  equations  of 
motion  for  the  bifocal  relay  mirror  spacecraft.  The  spacecraft  sketch  and  relative 
size  are  shown  in  Figure  1 .  It  consists  of  two  bodies,  transmit  telescope  and 
receive  telescope.  The  receive  telescope  rotates  with  respect  to  the  transmit 
telescope  about  a  single  axis.  Figure  2  depicts  a  typical  concept  of  operational 
employment  for  the  Bifocal  Relay  Mirror  Spacecraft.  The  receive  telescope 
acquires  and  tracks  the  source  laser,  while  the  transmit  telescope  tracks  the 
target.  The  motion  of  the  fast  steering  mirrors  accomplished  internal  alignment 
and  fine  pointing  of  the  beam  on  the  target. 

A.  COORDINATE  SYSTEM  DEFINITION 

The  coordinate  frame  descriptions  are  with  respect  to  Figure  3.  The 
center  of  mass  of  the  receive  telescope  is  on  the  rotation  axis.  Therefore,  the 
center  of  mass  of  the  system  is  fixed  during  the  relative  motion  of  the  receive 
telescope. 


Figure  1  Bifocal  Relay  Mirror  Spacecraft  sketch  (left)  and  relative  size  (right) 
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Figure  2  Concept  of  Operational  Employment  for  Bifocal  Relay  Mirror  Spacecraft 
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The  rigid  body  equations  of  motion  for  the  transmitter  and  receiver  are 
with  respect  to  the  inertial  frame.  Additionally,  the  dynamic  coupling  effects  from 
the  fast  steering  mirrors  are  neglected  due  to  the  small  mass  relative  to  the  entire 
structure. 

The  coordinate  system  fv  f2,  f3  is  fixed  in  the  transmit  telescope  with  £, 
axis  parallel  to  rotation  axis  of  the  receive  telescope,  z3  as  telescope  axis  and  z2 
is  normal  to  f,  and  z3  such  that  the  fv  f2,  z3  coordinate  system  is  right-handed 

and  mutually  orthogonal  frame.  The  origin  of  the  coordinate  system  is  at  center 
of  mass  of  the  transmit  telescope. 

The  coordinate  system  pv  p2,  p3  is  fixed  in  the  receive  system  with  pA 
axis  as  the  rotation  axis,  parallel  to  r,  axis,  p3  as  telescope  and  p2  is  normal  of 
,0,  and  p3  such  that  pvp2,p3  coordinate  system  is  right-handed  mutually 
orthogonal  frame  with  origin  at  the  center  of  mass  of  the  receive  telescope.  The 
coordinate  system  pv  p2,  p3  is  obtained  from  the  coordinate  system  zvz2,z3  by 

rotation  a  about  x-axis. 

The  relative  motion  of  the  transmit  fast  steering  mirror  with  respect  to 
transmit  telescope  is  represented  by  [iTx  rotation  about  £,  axis  and  (5Ty  rotation  f2 
axis.  The  relative  motion  of  the  receive  fast  steering  mirror  with  respect  to 
receive  telescope  is  represented  by  /?f  rotation  about  y a,  axis  and  rotation 

about  p2  axis.  The  fast  steering  mirror  rotations  are  assumed  to  be  small.  The 
equations  of  motion  of  the  system  are  written  in  the  coordinate  frame  rv  z2,  f3 . 


B.  EQUATIONS  OF  MOTION 

In  a  general  case,  rotational  equations  of  motion  for  a  given  body  about  an 
arbitrary  point  P  are  given  by 
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Equation  1:  Moment  and  Momentum  Balance  for  Rigid  Body 

Mp  =  fip-A  x  m  rc 

where 

Mp  =  total  sum  of  external  forces  about  P 

Hp  =  angular  momentum  of  the  body  about  P 

pc  =  total  vector  from  P  to  center  of  mass  of  the  system 

m  =  total  mass  of  the  body 

rc  =  velocity  of  center  of  mass  of  the  body 

If  the  point  P  is  the  center  of  mass,  then/?c=0,  then/W  =  H|w,  with  N 

being  used  to  indicate  the  inertial  reference  frame.  The  angular  momentum  of 
the  system,  Hs  can  be  written  as  follows: 

Equation  2:  System  Momentum  Balance 

HS=H  +  Hrel+Hw 
where 

Hs  =  total  angular  momentum  of  the  system. 

H  =  total  angular  momentum  by  neglecting  the  contribution  by  the 
relative  motion  of  the  receive  telescope  and  reaction  wheels  with  respect  to 
transmit  telescope 

Hrel  =  angular  momentum  due  to  relative  motion  of  the  receive  telescope 
Hw  =  angular  momentum  due  to  relative  motion  of  the  reaction  wheels. 

To  continue  the  momentum  derivation,  we  must  determine  the  inertia 
properties  of  the  system. 
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1.  Inertia 

Let  lT  be  the  inertia  matrix  of  the  transmit  telescope  about  its  center  of 


mass  in  coordinate  frame  fv  f2,  f3  ■ 

"  lT  -V  -V  " 

xx  xy  xz 

i  =  -r  r  -r 

t  yx  yy  yz 

-r  -r  r 

zx  zy  zz 

Let  the  vector  from  center  of  mass  of  the  transmit  telescope  to  the  center 
of  mass  of  the  spacecraft  is  given  by 

rBr  =xr?i+yrP2+zrP3 

Let  /',  be  the  inertia  matrix  of  the  receive  telescope  about  its  center  of 
mass  in  coordinate  frame  pv  p2,  p3 .  It  is  given  by 


f  p 

- lp ;  , 

x  y 

-r, 

x  z 

x'x' 

lp  , 

y  x 

lp,  , 

y  y 

-ip, 

y  z 

r, , 

z  X 

- ip :  . 

zy 

r, , 

z  z 

Next  we  transform  I'  to  f  f  f  frame.  The  transformation  matrix  from 

rv  f2,  f3  to  pv  p2,  p3  is  given  by 


P, 


CT 


1  0  0 
0  cos  (a)  sin(a) 
0  -sin(a)  cos(a) 


The  transformation  matrix  TCP  from  pv  p2,  p3  to  f  f  f  is  transpose 
Of  PCT .  /;  in  fv  f2,  f3  system,  lp ,  is  given  by  lp  =  TCp!'p  PCT 

Let  the  vector  from  center  of  mass  of  the  receive  telescope  to  the  center  of 
mass  of  the  system  is  given  by 
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The  inertia  matrix  of  the  spacecraft  about  its  center  of  mass  in  coordinate 
frame  fv  f2,  f3  is  given  by 


2  2 
y:  +z^ 

J  T  T 

-XrYr 

~XrZr 

r  2  2 

y  +z 

y  p  P 

-*PyP 

-*fyP " 

L  =1  +1  +m 

o  r  p  T 

-x  y 

tJ  t 

z2  +  x2 

T  T 

,  1 

+  m 

p 

X 

1 

z2  +  x2 

p  p 

l 

N 

-X  z 

T  T 

-y  z 

J  T  T 

x  +y 

T  J  T  J 

_  -XPZP 

1 

N 

2  2 
xn  +  y„ 

p  j  p  j 

2.  Angular  Velocities 

The  inertial  angular  velocity  of  the  transmit  telescope 
is N3r  =[®x  coy  ®z]  T .  The  relative  angular  velocity  of  the  receive  telescope 

with  respect  to  transmit  telescope  is  given  by Ta>p  =\a  0  0]  T .  The  inertial 

angular  velocity  of  the  receive  telescope  Ncop  =  NcoT  +  Ta>p  =  [<»x  +a  coy  «z]  1 . 

3.  Angular  Momentum 

The  terms  of  Equation  2  are  further  defined  below.  The  angular 
momentum  without  receive  telescope  relative  motion  is  given  by H  =  Isncdt  .  The 

relative  angular  momentum  is  given  by  Hrel  =  IP  rvp  ■ 

Through  substitution,  we  have 

Equation  3:  Refined  System  Momentum  Balance 

Hs=lsNtf+l/d>p+Hw 

4.  Kinematical  Relationships 

Using  Equation  1 ,  the  equation  of  motion  of  the  spacecraft  is  given  by 

Equation  4:  Rigid  Body  Moment  &  Momentum  Defined 

N  d  -  T  d  - 
M=  — H,  =  + 

dt  s  dt  s  s 

where 


N 


— Hs  =  rate  of  change  in  inertial  frame 


8 


-^-Hs  =  rate  of  change  in  transmit  telescope  frame 


It  should  be  noted  that  since  ls  is  function  of  a,  it  is  therefore  time 
dependent. 

a.  Transmitter  Kinematics 

The  momentum  balance  of  the  transmitter  is  as  given 
as:  nHt  =  lr  NooT  +  hw .  The  total  moment  applied  to  the  transmitter  is  given  by: 


N  d  -  N  d 

Mr  =  Mc  +  MG  =  -NH<=  -(lTN^+hw) 

•••=  ^{lT  ng>t  +  hw)+Nd>Tx  (/r  ng)t  +  hw ) 

-  =  (lTNo}T  +  fiw)+N^x(lrN^  +  hw) 

The  subscripts  C  and  G  represent  Control  and  Gimbal, 
respectively.  The  gimbal  torque  is  due  to  the  relative  motion  of  the  two 
telescopes.  The  receiver  momentum  is  given 

by NHP  =  I pNa>p  =  rCpl/CT{N3T  +  T3P)=  rC/7/Cr(/V  +5).  The  moment  applied 
to  the  receiver  is  given  by: 

N  d  -  N  d 

M  =-Mr  =  —nHp=  —{tCpI  pCT{N3T+d)) 

p  G  dt  dr  p 

■■■=  —{rCplpCT{Ncdr+d))+NarxTCpl  pCr{Nd)T+d ) 

■■■  =  (  TCplp  pCr  +  TCplp  PCT  )(NcoT+d)+  rCplp  pCT{Nd)T  +d) 

•••+'Vx  TCPI/CT{N3T  +cc) 

System-level  open-loop  control  input  is  achieved  through  system- 
level  dynamic  equilibrium,  as  follows. 
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Equation  5:  Open-Loop  Control  Torque 

Mc  =  -Mg  +{Mc+Mg)  =  Mp  +  Mr 

•••  =  {( TCplp  PCT  +  rCplp  pCr)(Ndr  +  a)  +  TCplp  pCr(Ndr  +d) 

■■■+Ncdrx  rcplp  pCT{NcbT  +  3))  +  {/r  Ndr  +  n3t  X  (/r  n3t  +  hw )} 

•  •  •  =  ip(N3T  +a)  +  !sNa>T  +  Ndr  X (/s ncot  +hw  +  TCplp pCTd)+  TCp\p PCT& 


In  later  chapters,  it  will  become  necessary  to  have  the  acceleration 
of  the  telescopes  in  order  to  determine  the  effects  of  flexibility.  Given  that 
attitude  and  angular  rates  are  observable  with  most  sensors,  the  above  equation 
may  be  solved  for  NcoT . 


Equation  6:  Angular  Acceleration  of  Transmitter 


N  r 


®t=ItotaL-\Mc  -  "cb'x{lTOTAL"cb'+hw  +  TCPI  pCTa )-!  (Naz  +a)~  rCpl/Cra) 


N  ~*T  . 


N  -*  T 


N  ~*T 


b.  Receiver  Kinematics 

The  inertial  angular  acceleration  of  the  receiver  is 

simply  No>p  =  NcoT  +a  . 


Equation  7:  Angular  Acceleration  of  Receiver 

n5)p  =  n5)t  +  5 

■■■  =  'total  \MC  -  Nco<  X  (lTOTAL  +  hw  +  TCpip  pC^)~  itN3'  +3)-  TCplp  pC^)  +  3 


5.  Wheel  Control  Laws 

Reaction  wheels  are  typically  controlled  by  a  proportional-derivative  (PD) 
controller  of  the  formH^  =2  kjqiEq4E  +  kida>iE .  This  control  law  is  developed  in 
Wie.  Where  qiE  is  the  ith  element  of  the  error  quaternion  and  cojE  is  the  error 

between  commanded  and  measured  angular  rate  about  the  ith  axis.  The  PD 
controller  used  in  this  work  is  of  the  same  form,  but  with  the  parameters  defined 
differently.  The  feedback  control  law  used  in  this  work  is  shown  in  Equation  8. 

Equation  8:  Proportional  Derivative  Controller 

HWi  =  kj(qiEq4E  +  t j  cOjE ) 
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Where  kj  is  the  controller  gain  and  r(  is  the  controller  time  constant.  This 
creates  a  dynamic  situation  of  the  form  Ip  +  k;TpiE  +  kjqiEq4E  =  0  .  The  bandwidth 
of  the  controller  goes  as  the  square  root  of  kj .  Hence,  increasing  the  bandwidth 

by  a  factor  of  two  requires  increasing  the  gain  by  a  factor  of  four. 

6.  Commanded  Trajectory 

The  commanded  trajectory  for  this  control  system  is  based  upon  initial 
conditions  for  the  joint  angle  motion  as  well  as  the  transmitter  acceleration  and 
rates.  The  initial  quaternion  for  the  transmit  telescope  is  given  by  qr(t0).  The 

initial  transmit  telescope  inertial  angular  velocity  is  given  by  NcoT(t0) .  The  initial 
joint  angle  is  given  by  a{t0 ) .  The  initial  joint  angle  is  given  by  a{t0 ) . 

Equation  9:  Initial  Conditions  for  Simulations 

0 
0 
0 
1 

"0 

N3T(t0)=  0  rad  I  sec 
0 

a{t0 )  =  0  rad 
a{t0 )  =  0  rad  /sec 

The  commanded  attitude  profile  is  based  upon  a  typical  engagement 
scenario.  The  inertial  acceleration  for  the  transmit  telescope  and  the  relative 
acceleration  for  the  receive  telescope  is  taken  from  the  existing  simulation.  The 
values  are  described  in  Equation  10. 

Equation  10:  Commanded  Acceleration  Profile  for  Simulations 

0.308sin(0.0079f) " 

Ntir(t)=  7.896sin(0.0126f)  xIO  5  rad/sec2 
-1 .097  sin(0.01 05f) 

d{t)  =  1.1 484 sin(0.0209f)x1 0  4  rad / sec2 
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III.  MULTI-BODY  FLEXIBLE  EQUATIONS  OF  MOTION 


A.  FLEXIBLE  EQUATIONS  OF  MOTION  -  GENERAL  CASE 

Given:  flexible  body,  with  a  cantilevered  attachment  at  point  A,  the 
equations  of  motion  may  be  stated  as  follows: 

Equation  11:  Cantilevered  Dynamics 


Where 

I  =  Internal  Nodes 
A  =  Attachment  Point 


M  &  K  =  Mass  &  Stiffness  Matrices  from  cantilevered  condition  in  finite 
element  analysis  program 


Equation  12:  Transfer  Nodes  to  Modal  Coordinates 


Where  qA  are  modal  coordinates  at  the  attachment  and  is  normalized 

for  unity  modal  masses,  i.e.c 1>etMiiOe=I.  The  terms  modal  coordinates  and 

generalized  coordinates  are  interchangeable  for  the  purposes  used  in  this 
document. 


Substitute  Equation  12  into  Equation  11  and  premultiply  by 


get: 


Equation  13:  Dynamics  Expressed  in  Modal  and  Attachment 
Coordinates 

I  Mae  \Qa\+  O)  0  j  )  _  f  /"/  j 

maet  mar\\ua\  |_o  oJKj  [Fj 

Where, 
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mae=oetm„or  +  oetmia 

Mar  =  ®rTM„®r  +  maior  +  &  rtmia  +  maa 

CO2  =  <£>/!<!!<&£ 

F,  =  ®eTF, 

Fa  =®rTF,+Fa 

Equation  13  is  general  in  nature,  and  specific  to  the  coordinates  of  the 
flexible  body.  The  next  series  of  equation  transforms  Equation  13  from 
coordinates  at  the  cantilevered  appendage  to  a  common  coordinate  system  with 
a  given  location  and  orientation.  The  selected  coordinate  system  is  defined  as 
the  axes  defined  by  the  transmit  telescope  coordinate  system,  located  at  the 
system  center  of  mass,  B*.  This  basis  is  called  {B},  as  it  is  the  body  coordinate 
system  use.  The  transmitter  basis  and  the  body  basis  differ  in  origin,  but  not 
orientation.  The  radius  vectors  from  the  B*  to  the  attachment  points  of  the 
transmitter  and  receiver  are  given  by: 

Equation  14:  Center  of  Mass  Locations 

rBV  =  xrP,  +  yrP2  +  zrP3 

rB’f=xM+yA+zA 

The  attachment  coordinates  are  expanded  to  divide  between  transmitter 
and  receiver  coordinates,  superscripted  by  %  and  p  respectively.  Translation  and 
rotation  are  indicated  by  lower  case  t  and  r  respectively.  For  each  flexible 

body {L/^ }  =  \UA  UrA^ .  Hence,  the  translational  and  rotational  attachment 
coordinates  may  be  expressed  in  terms  of  the  body  coordinates  as  follows. 

Equation  15:  Coordinate  Transformations 

ua\JCa  °P  CaRVu'} 

urA J  L  0  cJLo  I \\url  Lo  CA\\ur\ 

where 

" 0  zA  -yA~ 

[R]=  -zA  0 

_yA  -xA  0  _ 
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Substitute  Equation  15  into  Equation  13  and  premultiply  by 


7  0 
0  V 


to 


obtain 


Equation  16:  EOMs  in  Modal  and  Body  Coordinates 


M 

M 


AEO 


ARO . 


/ 

MaeoT 

where 


M AEO  -  MAE~T 
^ ARO  =  T  MArT 


'  Al  -  '  I 

F  =  TTF 

'  AO  11  A 


Qa 

0 


+ 


6)Z  0 

0  0 


u  = 


u 1 
uF 


\Qa 

u 


Al 


AO  _ 


B.  FLEXIBLE  EQUATIONS  OF  MOTION  -  TRANSMIT  TELESCOPE 

Flexible  Body  Equations  of  Motion  for  Transmit  Telescope  (Note:  A 
replaced  by  x) 


Coordinate  Transfer: 

{T}='CB{B} 


r 

_Cr 

1  2 

>  =  < 

lb  1 
|U2  1 

r 

00 

V _ 

CO 

-Cl 

_ J 

tCb  =  I 
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Equation  17:  Transformation  and  Mass  Matrices,  Transmit  Telescope 

"1  0  0  0  zr  -y~ 

0  1  0  -zr  0  xr 

r-^n  0  0  1  y  —x  0  [/  R~ 

L  rJ  0  0  0  1  0  0  |o  I 

0  0  0  0  1  0 

0  0  0  0  0  1 

~mr  0  0  0  0  0" 

0  mr  0  0  0  0 

,M  1  0  o  m,  0  0  0  JM,  0- 

1  ,Ri  0  0  0  /;,  /;  4  [o  /'. 

o  o  o  r  it  r 

xy  yy  yz 

o  o  o  r  r  r 

xz  yz  zz 


Mr  MrRr 
r]mt  RTzMTRT +r 


C.  FLEXIBLE  EQUATIONS  OF  MOTION  -  RECEIVE  TELESCOPE 

Flexible  Body  Equations  of  Motion  for  Receive  Telescope  (Note:  A 
replaced  by  p) 

Coordinate  Transfer: 

{R}=  PCB{B} 

/?.,  i  r  i  o  o  Vb, 

p2\=  0  cos(a)  sin(a)  <b2 
p3 \  0  -sin(a)  cos(a)  \b3 

"10  0 
PCB  =  0  cos(a)  sin(a) 

0  -sin(a)  cos  (a) 
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Equation  18:  Transformation  and  Mass  Matrices,  Receive  Telescope 


" PCB 

pcbr~ 

0 

PCB 

1 

0 

0 

0 

zp 

-yP 

0  1 

cos(a) 

sin(« 

) 

yps\n(a)-zpcos{a) 

~xp  sin(a) 

xp cos  (a) 

0  - 

-sin(a) 

cos  {a) 

ypcos(a)  +  zps\n(a) 

-x^cos  (a) 

~xps\n(a) 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

cos(a) 

sin(a) 

0 

0 

0 

0 

-sin(a) 

cos(a) 

0 

0 

0 

0 

0 


0 

mp 

0 

0 

0 

0 


0 

0 

m 

P 

0 

0 

0 


0 

0 

0 


r 

xx 

r 

xy 

ip 

xz 


o 

o 

0 


r 

xy 

r 

yy 

r 

yz 


0 

0 

0 


r 

xz 

£ 

£ 


M 


0 

r 


[mpRO]=ttpm 


pfJ" p 


'BT 


M  PCE 

p 


RTPPC 


BT 


Mp  pCb 


'BT 


M  PC  R 

p  p 


RTppc 


BT 


Mp  PC  Rp  + 


CBTlppCb 


Now,  we  have  Maro  for  transmit  and  receive,  we  need  Maeo  for  transmit 
and  receive.  See  Equation  16.  These  rigid-elastic  coupling  matrices  are 
available  from  most  of  the  widely  used  finite  element  analysis  programs,  such  as 
NASTRAN  or  l-DEAS. 
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IV.  STRUCTURE-CONTROLS  INTERACTION 


In  order  to  control  the  motion  of  flexible  structures  consisting  of  rigid  body 
and  elastic  motion,  the  coupling  between  rigid  body  motion  and  elastic  motion  is 
important.  Controlling  the  rigid  body  attitude  is  one  of  the  primary  objectives. 
Initially,  the  objective  of  achieving  fine  pointing  may  drive  the  control  design  to  a 
simple,  high  gain  controller.  That  is,  retain  the  commonly  used  PD  control 
architecture  with  higher  proportional  gain.  This  provides  more  control  authority 
based  to  control  attitude  errors.  If  the  designer  increases  the  proportional  gain 
by  a  factor  of  four,  the  control  bandwidth  doubles.  This  was  discussed  when 
Equation  8  was  presented.  If  this  high  gain  controller  is  within  the  structural 
frequency,  the  structural  modes  will  get  excited.  Given  the  rigid-body  model,  the 
gain  would  be  set  in  order  to  achieve  the  design-specific  pointing  accuracy. 
Ultimately,  that  is  how  the  PD  gains  would  be  adjusted  to  achieve  high  pointing 
accuracy.  The  question  is:  how  do  we  modify  the  controller  to  behave  favorably 
around  the  structural  frequencies? 

A.  COMPENSATOR  DESIGN 

Compensators  are  applied  after  the  rigid  body  control  is  applied.  The 
purpose  of  the  compensator  is  to  add  stability  to  the  system.  The  stability  is  with 
respect  to  gain,  phase,  or  a  combination.  Classical  gain-phase  stabilization 
techniques  are  of  primary  use  to  single-input  single-output  (SISO)  systems.  A 
rigid  or  flexible  body  is  said  to  be  gain-phase  stabilized  if  it  is  closed  loop  stable 
with  finite  gain  and  phase  margins.  These  classical  techniques  will  be  applied  to 
this  multi-input  multi-output  (MIMO)  system  [Ref  Wie], 

1.  Gain  Stabilization 

Gain  stabilization  of  a  flexible  body  attenuates  the  gain  at  desired 
frequencies,  to  provide  closed  loop  stability  regardless  of  phase  uncertainty.  A 
gain  stabilized  mode  has  a  finite  gain  margin,  but  is  closed  loop  stable  regardless 
of  the  phase  uncertainty. 
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2.  Phase  Stabilization 

Phase  stabilization  of  a  flexible  body  provides  phase  characteristics  at  the 
desired  frequency  to  obtain  closed  loop  damping  higher  than  the  passive 
damping  of  the  mode.  A  phase  stabilized  mode  has  a  finite  phase  margin,  but  is 
closed  loop  stable  regardless  of  the  loop  gain  uncertainty.  Since  the  dynamic 
model  has  no  passive  damping,  phase  stabilization  is  required  for  this 
application. 

B.  COMMON  CLASSICAL  COMPENSATOR  DESIGN 

Wie  presents  a  comprehensive  discussion  compensator  design  in  his 
book.  One  fundamental  concept  is  that  gain  stabilization  may  only  be  achieved  if 
the  corresponding  mode  has  some  passive  damping.  Since  no  modal  damping 
is  modeled,  phase  stabilization  techniques  must  be  employed.  The  following 
sections  highlight  some  common  compensator  designs. 

Control  theory  has  made  significant  advances  in  the  past  several  decades. 
Application  of  these  theories  has  been  difficult,  from  a  practicality  perspective. 
More  often  than  not,  classical  techniques  are  applied.  Performance  is 
“optimized”  through  trial  and  error  synthesis  using  frequency  domain  methods  or 
actual  performance  in  situ.  Due  to  the  complexity  of  the  dynamic  model, 
performance  will  be  compared  based  on  the  performance  of  the  plant  under 
various  compensation  techniques.  For  lower  performance  requirements,  simple 
first  order  filters  are  sufficient.  The  effect  of  the  first  order  filter  is  to  change  the 
controller  bandwidth.  Typically,  low  pass  filters  are  used  for  this  type  of  an 
application,  since  they  will  eliminate  the  high  frequency  content  of  the  controller, 
hence  not  exciting  the  structural  modes.  However,  for  higher  performance 
requirements,  such  as  placing  the  control  bandwidth  close  to  or  within  the 
structural  modes,  higher  order  filters  become  necessary. 

Wie  presents  the  generalized  second  order  filter,  using  the  s-domain.  He 
presents  it  as  a  natural  extension  of  the  classical  notch  or  phase  lead/lag  filter, 
and  is  based  on  pole-zero  patterns.  The  following  sections  employ  Wie’s 
descriptions  for  various  compensation  methods.  The  filter  is  of  the  form: 
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Equation  19:  Generalized  Second-Order  Filter,  Ref  Wie 


H(S)  = 


V-^+i 

a a 


-  +  - 


2Cds 


+  1 


a>„ 


(Or, 


P  P 

1.  Minimum-Phase  Lead  Filter 

A  phase-lead  filter  utilizes  Equation  19,  and  assigns^  =^z  =^p  >0  .  The 
maximum  phase  lead,  obtained  at coc  =  ^co/o~p  ,  is  determined  by: 

Equation  20:  Maximum  Phase  from  a  Phase  Lead  Filter,  Ref  Wie 

(2CV*V^z)  -(^p/^z-1)2 


6  =  cos  1 

r  max 


(2Cc  J  cop  I  coz)  +{(Opl(Oz-i)2 


,  dt  (Oc 


(Oz(Op 


The  gain  increases  or  decreases  at  high  frequencies,  and 

becomes  =40log10(«p/<yz)  dB  .  For  small  when  filter  poles  and  zeros  are 

close  to  the  imaginary  axis,  the  effective  lag  region  lies  between  coz  and  cop,  and 
the  maximum  phase  shift  approaches  ±180°.  In  practice,  ^pis  often  selected 
greater  than^z.  This  places  the  filter  poles  sufficiently  far  to  the  left  of  the 
imaginary  axis.  The  use  of  a  phase-lag  filter  with  a  large  cozlcop  ratio  greater 

than  two  should  be  avoided  from  a  practical  viewpoint.  The  phase  lag  filter  is 
shown  in  Figure  4. 
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Bode  Diagram 


Frequency  (rad/sec) 

Figure  4  Bode  Plot  of  Minimum-Phase  Lag  Filter 

2.  Minimum-Phase  Notch  Filter 

For  a  =a>z,  a  notch  filter  is  obtained.  The  minimum  gain  of  the  filter  is 
obtained  at  coc  =coz  =  cop  as Kx  =20log  10(^z/^p)  dB .  Both  phase  lead  and  lag 
occur  near coc.  For  the  notch  filter,  the  maximum  phase  lag  and  lead  occur  at  ca, 
and co2,  respectively,  where: 

Equation  21:  Location  of  Maximum  Phase  Lead  or  Lag  for  Notch 
Filter,  Ref  Wie 

^-  =  ^P+1-^(2C4+1)2-^ 

CLfc 

TT  =  + 1  +  ^(2(z(p  + 1)'  -1 

Wc 
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Because  —  and  —  depend  only  on  the  product  ££  ,  the  filter  damping 

0Jc 

ratios  determine  the  effective  notch  region.  Typical  notch  filter  damping  ratios 
are <^z  =  0,  and  £p  =  1 .  The  minimum  phase  notch  filter  is  shown  in  Figure  5. 


Bode  Diagram 


10'1  10°  101  102  103 

Frequency  (rad/sec) 

Figure  5  Bode  Plot  of  Minimum  Phase  Notch  Filter 

3.  Nonminimum-Phase  Notch  Filter 

Fo’ca>z=o)p,  and^p  >|^z|(^z  <0),  a  nonminimum-phase  notch  filter  is 

realized.  The  applicability  of  this  filter  is  when  passive  damping  does  not  exist 
(the  ideal  case).  The  conventional  notch  filter  cannot  be  used  in  this  case.  As 
the  pole-zero  combination  is  placed  further  from  the  imaginary  axis,  the 
robustness  of  the  filter  is  enhanced,  while  the  stability  margin  is  reduced.  The 
nonminimum  phase  notch  filter  is  shown  in  Figure  6. 
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Bode  Diagram 


10'1  10°  101  102  103 


Frequency  (rad/sec) 

Figure  6  Bode  Plot  of  Noninimum  Phase  Notch  Filter 

C.  CLASSICAL  COMPENSATOR  DESIGN 

Wie  presents  an  algorithm  for  designing  compensators  with  generalized 
filters.  The  basic  idea  is  to  synthesize  the  compensator  for  one  mode  at  a  time. 
The  design  algorithm  is  summarized  below: 

1.  Select  Control  Bandwidth 

The  driving  control  bandwidth  is  one  of  the  key  parameters  of  the  design. 
This  may  be  driven  by  performance,  noise  sensitivity,  control  authority,  etc. 

2.  Rigid  Body  Mode  Compensation 

Rigid  body  control  is  typically  achieved  through  a  proportional-derivative 
(PD)  feedback  controller.  The  format  of  the  control  law  is  given 
by:  u  =  -K(  1  +  Ts)y  =  ~(KP  +  KDs)y  . 
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3.  Flexible  Mode  Compensation 

The  types  and  degrees  of  the  closed  loop  behavior  depend  on  the  relative 
spectral  separation  of  the  modes.  Determine  if  simple  gain  stabilization  of  all  the 
flexible  modes  are  possible.  If  not,  determine  the  needed  lead  or  lag  angles  foe 
each  mode.  Synthesize  the  appropriate  structural  filter  for  each  destabilizing 
mode  one  by  one,  using  the  various  second-order  filters. 

4.  Design  Iteration 

Repeat  the  process  to  compromise  the  interactions  between  each 
compensator.  A  few  iterations  using  a  software  package,  or  analytical  model,  will 
result  in  a  more  streamlined  design  process. 
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V. 


DEVELOPMENT  AND  VALIDATION  ANALYTICAL  MODEL 


A.  FLEXIBILITY  SIMULATION  ARCHITECTURE 

Now  that  the  flexible  telescopes  are  modeled  dynamically,  a  method  of 
simulation  must  be  developed.  The  approach  used  is  the  feedback  torque  effect. 
This  method  is  presented  in  simple  terms  in  Figure  7.  In  essence,  the  rigid  body 
dynamics  are  unchanged.  However,  the  state  vector  determines  the  base 
excitation  (acceleration)  which  creates  a  feedback  torque.  This  feedback  torque 
is  combined  with  the  control  torque  to  affect  the  rigid  body  dynamics. 


T - O - . 


Flexible  Body 

Dynamics 

Figure  7  Flexibility  Model,  After  Ref  Agrawal 

In  Figure  7,  T  represents  the  command  input  (Torque),  and  0b  represents 
the  state  vector.  Using  Equation  16  as  a  reference,  the  general  equation  that  we 
are  solving  is  of  the  formAx  +  Bx  =  F ,  where  x  =  [q,  U]T.  The  most  practical  and 
mathematically  stable  method  of  solving  this  system  of  equations  is  to  determine 
the  rigid  body  state  and  work  backwards  to  determine  the  reaction  torque  which 
results  from  the  flexibility.  This  is  shown  in  the  next  series  of  equations.  If  we 
start  with  Equation  16  as  a  base  equation,  the  known  values  are  the  mass  and 
stiffness  matrices,  as  well  as  0  and  qA .  We  may  restate  Equation  16  as  two 
differential  algebraic  equations,  as  in  Equation  22. 


Rigid  Body 

Dynamics 

w 
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Equation  22:  Application  of  Equation  16  for  Determining  Simulation 
Algorithm 

lqA+MAEOU  +  a2qA=0 

MaEO  qA+  M ARoU  =  TA 

Through  reorganization,  Equation  22  may  be  restated  as 

Equation  23:  Determining  Reaction  Torque 

qA  =  -mAec>u  ~  ®  qA 

T reaction  =  AEO  ^A  +  ^ ARO^) 

At  first,  it  isn’t  intuitive  why  the  second  equation  has  the  negative  sign. 
Without  that  negative  sign,  the  system  is  unstable.  With  the  negative  sign,  the 
system  behaves  as  it  should.  In  plain  verbiage,  the  negative  sign  is  there 
because  it  is  a  reaction  torque,  the  behavior  of  which  opposes  that  of  an  applied 
torque.  The  generalized  coordinates  are  initialized  at  zero  when  the  simulation 
begins.  The  stable  algorithm  for  determining  the  reaction  torque  is  shown  in 
depicted  in  Figure  8. 


Figure  8  Determining  Feedback  Torque  from  Rigid  Body  Motion 

B.  DYNAMICAL  SIMULATION  DEVELOPMENT 

The  base  acceleration  is  determined  using  Equation  6  and  Equation  7  for 
the  transmitter  and  receiver,  respectively.  The  block  diagram  in  Figure  9  shows 
how  to  determine  the  inertial  acceleration  of  each  body  from  the  observed  states, 
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using  Equation  6  and  Equation  7.  The  block  diagram  in  Figure  10  determines  the 
feedback  torque  due  to  flexibility,  using  Equation  23. 


12  out 


Figure  9  SIMULINK  Model  for  Determining  Rigid  Body  Acceleration  of  Transmitter 

and  Receiver,  using  Equation  22  and  Equation  23 
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HI 

Figure  10  SIMULINK  Model  for  Determining  Torque  Due  to  Acceleration  of  Flexible 

Bodies,  using  Equation  23 

C.  VALIDATION  OF  DYNAMICAL  MODEL 
1.  Developing  Flexible  Bodies 

The  MSC  software  suites  of  PATRAN  and  NASTRAN  were  used  for 
modeling  the  geometry  and  determining  the  modes  of  each  telescope.  The 
objective  was  to  have  a  first  mode  for  the  transmitter  near  10  Hz,  with  the 
receiver  being  about  slightly  higher,  due  to  less  mass  and  the  same  stiffness. 
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This  requirement  was  changed  to  around  1  Hz  for  the  purposes  of  propagating 
an  answer.  Lower  modal  frequencies  take  less  time  to  propagate. 

a.  PATRAN  Telescope  Geometry 


Transmitter 

Receiver 

Inner  Diameter  (m) 

1.573 

1.573 

Outer  Diameter  (m) 

2.360 

2.360 

Length  (m) 

4.000 

4.000 

Mass  (kg) 

2267.6 

972.36 

Density  (kg/mA3) 

233.23 

100.01 

Young's  Modulus  (Pa) 

1.00E+08 

1.00E+08 

Poisson's  Ratio 

0.3 

0.3 

Figure  1 1  Flexible  Body  Parameters  from  Geometry  in  PATRAN 

b.  NASTRAN  Normal  Modes 

The  resultant  modes  from  NASTRAN  are  shown  in  Figure  12. 
During  simulation,  the  impacts  of  flexibility  were  impractical  for  such  high 
frequencies.  Therefore,  the  modes  were  modified  with  a  flexibility  factor.  It  was 
determined  that  reducing  the  modes  by  a  factor  of  ten  resulted  in  a  manageable 
simulation.  The  modes  used  for  simulation  are  shown  in  Figure  13. 


Mode 

Transmitter 

Receiver 

1 

8.59 

13.09 

2 

10.31 

15.72 

3 

16.73 

25.52 

4 

18.52 

28.24 

5 

27.64 

42.16 

6 

39.58 

60.35 

7 

47.99 

73.18 

8 

49.67 

75.75 

9 

51.92 

79.18 

10 

67.57 

103.04 

Figure  12  Receiver  and  Transmitter  Modes  from  NASTRAN  (Hz) 
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Mode 

Transmitter 

Receiver 

1 

0.86 

1.31 

2 

1.03 

1.57 

3 

1.67 

2.55 

4 

1.85 

2.82 

5 

2.76 

4.22 

6 

3.96 

6.04 

7 

4.80 

7.32 

8 

4.97 

7.58 

9 

5.19 

7.92 

10 

6.76 

10.30 

Figure  13  Receiver  and  Transmitter  Modes  Used  in  Simulation  (Hz) 

c.  Modeling  Rigid-Elastic  Coupling 

The  rigid-elastic  coupling  matrix  was  not  taken  from  a  finite  element 
analysis  program.  Instead  was  chosen  arbitrarily,  from  visual  inspection  of  the 
displayed  modes  in  the  simulation  mode  of  PATRAN.  However,  the  method 
used  could  easily  be  substituted  by  NASTRAN  output.  The  arbitrary  method 
used  involved  selecting  a  unit  vector  for  each  mode  considered,  and  then 
multiplying  this  by  a  coupling  factor.  The  unit  vectors  represent  the  amount  of 
the  coupling  that  is  occurring  between  an  axis  of  rotation  and  a  generalized 
coordinate.  The  first  two  modes  are  bending  modes  about  the  x  and  y-axis 
respectively.  The  third  mode  is  a  twisting  mode  about  the  z-axis.  The  fourth  and 
fifth  modes  are  combinations  of  these  effects.  The  rigid-elastic  coupling  matrix 
was  adjusted  such  that  the  maximum  coupling  occurred  from  the  first  mode,  and 
the  minimum  coupling  occurred  from  the  last  mode  considered.  The  coupling 
factors  varied  from  a  maximum  of  10%  to  a  minimum  of  1%  of  the  maximum 
principal  moment  of  inertia,  and  are  shown  in  Figure  14. 


Mode 

Coupling  Factor 

1 

7000 

2 

5425 

3 

3850 

4 

2275 

5 

700 

Figure  14  Rigid-Elastic  Coupling  Factors  for  Simulation 
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The  rigid-elastic  coupling  matrix  is  presented  in  the  following,  where 

MaEO  =  M pEO  =  MtEO  ■ 

Maeo  =  CouplingFactorMatrix  *  CouplingVectors 


7000 

0 

0 

0 

0  " 

1 

0 

0 

0 

5425 

0 

0 

0 

0 

1 

0 

0 

0 

3850 

0 

0 

0 

0 

1 

0 

0 

0 

2275 

0 

0.7071 

0 

0.7071 

0 

0 

0 

0 

700 

0.7071 

0.7071 

0 

7000 

0 

0 

0 

5425 

0 

0 

0 

3850 

kg 

-m2 

1608.67  0  1608.67 

494.975  494.975  0 


2.  Observing  System  Free-Free  Behavior 

Nominally,  all  generalized  coordinates  are  assigned  an  initial  condition  of 
zero.  In  order  to  demonstrate  the  uncontrolled  dynamics,  generalized 
coordinates  are  given  random  initial  conditions,  and  no  actuation  is  applied.  The 
telescope  joint  angle  motion  and  spacecraft  attitude  actuators  (reaction  wheels) 
are  set  to  zero  input.  This  is  equivalent  to  stretching  and  compressing  various 
springs  in  an  undamped  multi-body  spring-mass  system  and  letting  them  go. 
The  purpose  is  to  observe  the  free-free  behavior  of  the  system. 

Figure  15  was  created  with  random  initial  conditions  for  the  generalized 
coordinates.  Additionally,  the  initial  angular  rate  of  the  transmitter  was  [0,0,0] 
rad/sec.  There  was  no  commanded  angular  acceleration  of  the  transmitter  or 
receiver,  and  there  was  no  initial  angular  momentum.  The  purpose  of  this  figure 
is  to  demonstrate  the  complex  frequency  content  of  the  response.  The  order  of 
magnitude  of  the  difference  in  angular  rates  is  not  to  be  compared  to  future 
graphs.  Since  the  initial  conditions  are  random,  any  such  comparison  would  be 
meaningless. 
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Difference  in  Angular  Rates 


x  10'6 


Figure  15  Uncommanded  Free-Free  Behavior,  Angular  Rates 

3.  Eigenvalues  of  the  Free-Free  System 

Fundamentally,  the  free-free  behavior  of  the  system  is  dependent  on  the 
joint  angle.  This  assumes  that  the  rate  of  motion  has  no  effects  on  the  free-free 
modal  behavior.  Since  the  input  has  no  effect  on  the  modes,  the  base  equation 
is  shown  in. 

Equation  24:  Free-Free  Dynamics,  No  External  Forces 


/ 

0 

M AEO 

Qr 

0)r2 

0 

o“ 

V 

0 

0 

/ 

M AEO 

< 

*P 

> 

0 

< 

0 

-< 

Ip 

>  =  < 

0 

Maeot 

MAEOT 

MARO_ 

u 

0 

0 

0 

u 

0 

When  this  is  put  into  the  eigenvalue  format,  we  have  Equation  25. 
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Equation  25:  Free-Free  Dynamics,  Eigenvalue  Solution 

det  {/l2/ -/VrXj  =  0 
where 


1 

0 

^  AEO 

M  = 

0 

/ 

M AEO 

AEO 

T 

M AEO  M aro 

COT2 

0 

o" 

K  = 

0 

O 

CM 

0 

0 

0 

A  simulation  was  created  to  determine  the  eigenvalues  of  the  system  at  a 
family  of  joint  anglesae[0  This  simulation  solves  Equation  25,  and  is 

shown  in  Figure  16. 


Figure  16  SIMULINK  Model  for  Determining  Eigenvalues  of  the  Free-Free  System, 

using  Equation  25 


The  eigenvalues  are  sensitive  to  the  coupling  factors  for  the  rigid-elastic 
coupling.  As  the  coupling  factors  are  reduced,  the  system  free-free  frequencies 
approach  the  modes  of  each  single  body.  These  free-free  modes  are  what  must 
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be  compensated  for  in  the  control  system.  For  the  coupling  factors  used,  the 
free-free  system  frequencies  are  shown  in  Figure  17.  The  column  identified  as 
“without  coupling”  uses  coupling  factors  of  zero.  In  all  columns,  the  rigid  body 
modes  are  omitted,  as  the  frequencies  are  zero  for  these  modes. 


Frequencies  (rad/sec) 

with  coupling 

without  coupling 

List  of  Single 
Body  Modes 

0 

5.3949 

5.3949 

0 

6.4756 

6.4756 

0 

8.2272 

8.2272 

6.947 

9.8752 

9.8752 

8.349 

10.5144 

10.5144 

11.374 

11.637 

11.637 

13.737 

16.0343 

16.0343 

17.299 

17.3696 

17.3696 

17.426 

17.7463 

17.7463 

26.411 

26.4884 

26.4884 

Figure  17  Free-Free  Modes  with  and  without  Rigid-Elastic  Coupling 

D.  SENSITIVITY  OF  PLANT  BEHAVIOR  TO  UNCOMPENSATED  PD 
CONTROLLER  BANDWIDTH 

This  section  shows  the  interaction  between  the  PD  controller  bandwidth 
and  the  flexible  plant,  without  any  compensator  applied.  This  is  accomplished  by 
turning  off  all  compensators  in  the  model,  holding  the  time  constant  unchanged, 
and  steadily  increasing  the  controller  gain. 

Fundamentally,  designing  a  control  law  for  a  system  with  flexible  bodies 
seeks  high  bandwidth  control  to  achieve  fine  pointing,  while  not  exciting  flexible 
modes.  When  the  original  rigid  body  gains  for  the  PD  controller  were  applied  to 
the  flexible  system,  they  were  too  low.  The  attitude  and  rate  errors  which 
occurred  were  due  to  a  low  controller  gain,  and  not  vibration.  To  achieve  higher 
pointing  accuracy,  the  controller  gains  were  increased  while  the  time  constants 
were  unchanged. 

The  original  rigid  body  controller  had  gains  around  5000  and  time  constant 
around  0.3  sec.  The  time  constant  was  adjusted  to  r  =  0.01  sec  in  order  to  lower 

the  damping  of  the  rigid  body  controller.  Then,  the  gains  were  adjusted  to  find 
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the  points  of  instability  for  various  modes.  The  highest  stable  gain  was 
determined,  as  well  as  the  gain  which  caused  instability  for  the  first  three  modes. 
The  adjusted  values  of  K  increase  the  bandwidth  of  the  controller,  and 
demonstrate  structure-controls  interaction. 

The  figures  that  follow  show  the  results  of  the  uncompensated  flexible 
body  for  increasing  gains.  There  are  three  plots  in  each  figure.  The  top  plot  is 
the  difference  in  quaternions.  This  has  no  physical  meaning,  but  does  allow  the 
user  to  see  the  rate  at  which  error  is  occurring.  The  middle  plot  is  the  difference 
in  angular  rates:  commanded  rate  minus  actual  rate.  The  bottom  plot  is  the 
pointing  error,  in  radians:  angle  subtended  between  the  actual  attitude  and  the 
commanded  attitude. 

In  Figure  21,  the  angular  rate  errors  are  on  the  order  of  1e-6  rad/sec,  and 
the  pointing  error  continues  to  grow.  This  is  an  effect  of  a  low  gain.  In  Figure  22, 
Figure  23  and  Figure  24,  the  angular  rate  errors  continue  to  decrease,  and  the 
onset  of  angular  error  is  delayed.  Finally,  in  Figure  25  there  is  no  angular  error 
throughout  the  simulation.  Now  the  gain  is  increased  until  the  plant  cannot  be 
controlled.  The  rate  error  continues  to  decrease  in  Figure  26  and  Figure  27. 
However,  in  Figure  28,  the  system  instability  is  obvious. 

Through  trial  and  error,  the  gains  at  which  each  axis  want  unstable  were 
determined.  Figure  29  shows  the  behavior  at  the  highest  stable  gain,  50.3e6.  At 
this  gain,  the  angular  rate  error  is  on  the  order  of  5e-10  rad/sec.  Uncompensated 
gains  above  this  level  result  in  rate  errors  on  the  order  of  5e-4  rad/sec,  one 
million  times  higher  than  the  highest  stable  gain.  The  first  instability  occurs  at 
50.4e6,  as  shown  in  Figure  30.  The  second  axis  is  unstable  at  85.7e6,  as  shown 
in  Figure  31 .  The  third  axis  is  unstable  at  93.6e6,  as  shown  in  Figure  32.  For  all 
of  these  simulations,  the  initial  pointing  error  was  10  nanoradians.  That  causes 
the  initial  oscillation. 
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Gain 

Max  Rate  Error 
at  20  seconds 
(rad/sec) 

Angle  Error  at 
20  seconds 
(rad) 

5,000.00 

5.E-06 

3.E-05 

20,000.00 

1.E-06 

8.E-06 

80,000.00 

3.E-07 

2.E-06 

320,000.00 

1.E-07 

3.E-07 

1,280,000.00 

2.E-08 

0 

5,120,000.00 

5.E-09 

0 

20,480,000.00 

1.E-09 

0 

50,300,000.00 

5.E-10 

0 

50,400,000.00 

5.E-04 

4.E-04 

93,600,000.00 

7.E-04 

6.E-04 

Figure  18  PD  Controller  Gain  and  Corresponding  Errors  for  Flexible  Spacecraft, 

Without  Compensator 


Gain 

Compensator 

Type 

Omega_z 

Omega_p 

Zeta_p 

Max  Rate  Error  at 
20  seconds 
(rad/sec) 

Angle  Error  at 
20  seconds 
(rad) 

51,000,000 

NMP  Notch 

first  three  modes 

-0.00001 

1 

5.00E-08 

0 

60,000,000 

NMP  Notch 

first  three  modes 

-0.00001 

0.7 

4.00E-07 

0 

55,000,000 

NMP  Notch 

first  three  modes 

-0.00001 

0.7 

3.00E-07 

0 

Figure  19  PD  Controller  Gain  and  Corresponding  Errors  for  Flexible  Spacecraft, 

Nonminimum  Phase  Notch  Filter 


Gain 

Compensator 

Type 

Omega_z 

Omega_p 

Zeta_p 

Max  Rate  Error  at 
20  seconds 
(rad/sec) 

Angle  Error  at 
20  seconds 
(rad) 

51,000,000 

Phase  Lag 

14 

6 

1 

1.5 

5.00E-10 

0 

60,000,000 

Phase  Lag 

14 

6 

1 

1.5 

4.00E-10 

0 

86,000,000 

Phase  Lag 

14 

6 

1 

1.5 

3.00E-10 

0 

95,000,000 

Phase  Lag 

14 

6 

1 

1.5 

3.00E-10 

0 

KBB 

Phase  Lag 

14 

6 

1 

1.5 

3.00E-10 

0 

luamm 

Phase  Lag 

14 

6 

1 

1.5 

2.00E-10 

0 

m&mm 

Phase  Lag 

14 

6 

1 

1.5 

2.00E-10 

0 

iwmm 

Phase  Lag 

14 

6 

1 

1.5 

1.80E-10 

0 

|  160,000,000| 

Phase  Lag 

14 

6 

1 

1.5 

1.50E-10 

0 

mm&m 

Phase  Lag 

14 

6 

1 

1.5 

1.00E-10 

0 

Phase  Lag 

20 

6 

0.3 

0.33 

1.20E-10 

0 

|  258,000,000| 

Phase  Lag 

20 

6 

0.3 

0.35 

1.20E-10 

0 

mmxm 

Phase  Lag 

20 

6 

0.2 

0.21 

1.20E-10 

0 

immm 

Phase  Lag 

20 

6 

0.2 

0.21 

5.00E-1 1 

0 

Figure  20  PD  Controller  Gain  and  Corresponding  Errors  for  Flexible  Spacecraft, 

Phase  Lag  Filter 
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Difference  in  Quaternions 


x  iq6  Difference  in  Angular  Rates 


-  x-axis 

-  y-axis 

-  z-axis 

Figure  21  Uncompensated  Behavior,  K  =  5000 


x  io"6  Difference  in  Quaternions 


x  10"6 


Angle  Error  (rad) 


-  x-axis 

-  y-axis 

-  z-axis 

Figure  22  Uncompensated  Behavior,  K  =  20,000 
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Difference  in  Quaternions 
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0  2 
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8  10  12 

14 

16 

18 

20 

x  10'7 

Difference  in  Angular  Rates 

x-axis 

y-axis 

z-axis 


Figure  23  Uncompensated  Behavior,  K  =  80,000 
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Figure  24  Uncompensated  Behavior,  K  =  320,000 
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Figure  25  Uncompensated  Behavior,  K  =  1,280,000 
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Figure  26 
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Figure  27  Uncompensated  Behavior,  K  =  20,480,000 
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Figure  28  Uncompensated  Behavior,  K  =  81 ,920,000 
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Figure  29  Uncompensated  Behavior,  K  =  50,300,000,  the  last  stable  gain 
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Figure  30  Uncompensated  Behavior,  K  =  50,400,000,  two  axes  unstable 
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Figure  31 


Figure  32 
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VI.  COMPENSATOR  APPLICATION  TO  ANALYTICAL  MODEL 


Based  on  the  compensators  described  in  Chapter  IV,  two  applicable 
designs  were  selected  for  application  to  this  control  system.  The  Nonminimum- 
Phase  Notch  Filter  and  the  Minimum-Phase  Lag  Filter  were  chosen  for  this 
undamped  system.  The  conventional  Minimum  Phase  Notch  Filter  is  not 
applicable  since  there  is  no  passive  system  damping  for  the  flexible  modes.  The 
structural  filters  used  in  industry  are  of  the  Minimum-Phase  Lead  variant.  This 
type  of  filter  was  attempted.  However,  no  performance  increases  were  realized. 
The  lowest  error  achieved  without  a  compensator  occurred  at  K  =  50.3e6.  At  this 
gain,  the  angle  error  is  numerically  “zero”.  Although,  it  is  not  zero,  it  is  below  the 
threshold  of  which  MATLAB  is  able  to  take  the  arccosine.  The  rate  error  at  the 
end  of  the  simulation  was  approximately  5e-10  rad/sec.  The  objective  of  the 
compensator  is  to  reduce  this  rate  error,  thereby  reducing  the  impact  of  structural 
vibrations  on  the  pointing  error  budget. 

A.  NONMINIMUM  PHASE  NOTCH  FILTER 

The  first  method  used  was  the  nonminimum  phase  notch  filter.  The 
compensators  were  designed  about  each  identified  free-free  frequency.  Initially, 
one  was  applied.  By  the  end,  three  compensators  were  applied,  with  the 
purpose  of  suppressing  the  first  three  modes.  The  best  performance  achieved 
was  on  the  order  of  3e-8  rad/sec.  Unfortunately,  this  is  worse  than  the 
uncompensated  system  before  the  onset  of  instability.  The  results  are  shown 
below. 
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Figure  33  Notch  Filter,  K  =  51  e6,  <^=-0.00001,^=1.0 
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Figure  34  Notch  Filter,  K  =  60e6,  <^z  =  -0.00001,  =  0.7 
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Figure  35  Notch  Filter,  K  =  55e6,  <^z  =  -0.00001,  =  0.7 

B.  MINIMUM  PHASE  LAG  FILTER 

The  second  method  used  was  the  minimum  phase  notch  filter.  The 
compensators  were  designed  about  the  identified  free-free  frequency.  The  best 
performance  achieved  was  on  the  order  of  5e-1 1  rad/sec.  This  is  a  factor  of  ten 
reduction  in  angular  rate  error. 

Initially,  a  highly  damped  filter  was  used.  With  <^z=1  and  ^p=1 .5  ,  the  gain 

was  steadily  increased  from  Figure  36  through  Figure  45.  The  rate  error  in 
Figure  45  was  reduced  to  1e-10  rad/sec.  At  this  point,  with  a  gain  of  258e6,  the 
error  has  been  reduced  by  a  factor  of  five. 

After  this,  the  filter  was  lightly  damped.  The  zeros  had  damping  ratios 
between  0.2  and  0.3.  The  poles  had  damping  ratios  1 .05  times  that  of  the  zeros. 
Different  cases  are  presented  in  Figure  46  through  Figure  49.  The  rate  error  in 
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Figure  49  was  reduced  to  5e-1 1  rad/sec.  At  this  point,  with  a  gain  of  500e6,  the 
error  has  been  reduced  by  a  factor  of  ten. 
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Figure  36  Phase  Lag  Filter,  K=51e6,  <^z=1,  <TP=1-5,  coz= 14,  cop= 6 
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Figure  39 


Figure  40 
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Phase  Lag  Filter,  K=95e6,  <^z=1,  ^=1-5,  ®z=14,  cop= 6 
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Figure  41 


Figure  42 
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Phase  Lag  Filter,  K=115e6,  <^=1,  £p=1-5,  ryz=14,  <yp=6 
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Figure  43 
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Figure  44  Phase  Lag  Filter,  K=160e6,  <^z=1,  <^p=1.5,  coz= 14,  cop= 6 
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Figure  45 


Figure  46 
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Figure  47  Phase  Lag  Filter,  K=258e6,  <^=0.33,  ^=0.35,  az=20,  cop=6 
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Figure  48  Phase  Lag  Filter,  K=258e6,  Cz=0.2,  <^p=0.21,  oz=20,  cop=  6 
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Figure  49  Phase  Lag  Filter,  K=500e6,  £,=0.2,  £P=0.21,  coz= 20,  cop=  6 
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VII.  CONCLUSION 


The  design  requirements  for  the  Bifocal  Relay  Mirror  spacecraft  include 
controlling  jitter  at  the  nanoradian  level.  Typically,  tight  pointing  requirements 
require  high  structural  stiffness,  at  the  cost  of  increasing  the  on-orbit  mass.  To 
accomplish  this,  while  minimizing  the  mass  of  the  spacecraft,  the  structure  will 
have  some  inherent  flexibility.  These  flexible  modes  will  interact  with  the  pointing 
control,  hence  affecting  the  payload  performance.  The  compensator  design 
conducted  in  this  thesis  achieved  order  of  magnitude  improvements  in  controlling 
the  rate  error,  hence  jitter.  This  thesis  presented  the  rigid  body  and  flexible  body 
dynamics,  and  compares  the  uncompensated  and  compensated  performance  of 
the  spacecraft  using  a  flexible  plant.  Classical  compensator  designs  were 
reviewed  and  applied. 

The  plant  in  this  work  was  nonlinear.  No  attempts  were  made  to  linearize 
the  plant  about  a  set  of  conditions,  or  for  a  given  set  of  command  inputs.  The 
instabilities  were  found  completely  by  trial  and  error.  One  alternate  approach  for 
determining  the  destabilizing  gain  would  be  to  command  a  rotation  about  one  of 
the  principal  axis.  Even  with  a  time  variant  parameter,  this  would  provide  a  linear 
plant.  This  would  allow  the  control  designer  to  use  classical  controller  design, 
such  as  root  locus  techniques,  to  determine  the  unstable  gains.  This  thesis  did 
no  such  linearization,  as  a  primary  task  was  to  establish  a  baseline  nonlinear 
simulation  to  which  others  could  be  compared. 

This  thesis  showed  that  classical,  high  gain  PD  controllers  are  appropriate 
for  rigid  body  control.  However,  when  the  controlled  plant  includes  flexible 
modes,  uncompensated  high  gain  PD  control  causes  extreme  deviation  due  to 
unmodeled  dynamics.  The  compensators  used  in  this  formulation  validate  their 
application  on  flexible  spacecraft  with  high  pointing  accuracy,  such  as  the  Bifocal 
Relay  Mirror.  The  relative  motion  of  the  two  telescopes  is  slow  in  this  case.  No 
effort  was  made  to  adjust  the  compensator  based  on  the  joint  angle  motion 
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characteristics.  Classical  compensator  synthesis  was  a  sufficiently  accurate 
method  for  controlling  the  feedback  torque  due  to  flexibility. 

Through  classical  second-order  compensators,  the  angular  rate  error  was 
decreased  by  a  factor  of  ten.  Nonminimum  phase  notch  filters  and  phase  lag 
filters  were  used.  Ultimately,  the  phase  lag  filters  provided  the  best  performance. 
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